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The phase diagram of the half-filled spinless Holstein model for elec- 
trons interacting with quantum phonons is derived in three dimensions ex- 
tending at finite temperature T a variational approach introduced for the 
one-dimensional T = case. Employing the variational scheme, the spectral 
and optical properties of the system are evaluated in the different regimes 
that characterize the normal and ordered state. The effects of the charge- 
ordering (CO) induce a transfer of spectral weight from low to high energies 
in the conductivity spectra, as the temperature decreases or the strength of 
the electron-phonon (el — ph) interaction increases. The inclusion of effects of 
lattice fluctuations is able to smooth the inverse square-root singularity ex- 
pected for the case of the mean-field approach and determines a subgap tail 
absorption. Moreover, in the weak to intermediate coupling regime, a two- 
component structure is obtained within the CO phase at low frequency: the 
remnant Drude-like term and the incipient absorption band centered around 
the gap energy. 
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In the last years there has been a renewed interest in charge density wave (CDW) 
materials. 1 The transition to a CO phase is common to a wide range of compounds, 2 includ- 
ing quasi-one- dimensional organic conductors, 1 ' 3 dichalcogenides, 4,5 molybdenum bronzes 2,6 
and A-15 materials. 7 Moreover, the CO transition has been associated with the stripe density 
wave order in some cuprates 8 and nickelates, 9 ' 10 with the commensurate or incommensurate 
charge ordering in manganites. 11 

The ordered phase generally evolves out of a metallic phase. In the weak el — ph coupling 
regime the transition is the well-understood CDW instability of the Fermi liquid. 1 The 
related equations predict a ratio of T = gap A to the ordering temperature Tco close 
to the BCS value A/ksT C o = 3.52. However in quasi-one-dimensional materials this ratio 
assumes larger values because of critical fluctuations in low dimensionality. Furthermore 
values of the ratio larger than 10 have been measured in quasi-two-dimensional and three- 
dimensional materials. 4 ' 5,7 ' 12-14 In Ni compounds 12 A/ksTco — 13 and in manganites 14 a 
ratio of order of 30 have been deduced by the energy gap fitted by optical measurements. 
This suggests that a strong coupling CO phase develops in novel materials and that the 
lattice degrees of freedom can be important for stabilizing the ordered state. 15 A very useful 
approach to investigate the properties of the ordered state is to study the optical absorption 
of CO materials. 1,15-18 In particular the lattice fluctuation effects 19-21 and the response of 
the system in the strong coupling CO regime 12 ' 15 ' 22 have been studied. Indeed the inverse 
square-root singularity expected for the case of a static distorted lattice can be removed by 
the fluctuation effects giving rise to a subgap tail absorption. 20 Furthermore the presence of a 
polaronic peak in the MIR band of Ni compounds observed above and below the transition 
temperature can be explained only within a theory valid for intermediate to strong coupling 

19 1^ 

regime. ' 

Theoretically one of the simplest models of electrons interacting with phonons and de- 
veloping a CO state is the Holstein model. 23 The CDW correlations and the transition 
temperature have been studied within this model by Monte Carlo simulations 24-28 showing 
that the ordered phase is stable at half-filling. Recently, using the dynamical mean-field 
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theory (DM FT), the issue of the non-perturbative coupling scenario of the CO phase and 
of the origin of the large ratio of the zero temperature gap to the transition temperature has 
been addressed. 29-31 The results point out that the scattering of electrons by the phonons 
involved in the lattice distortions represents the crucial effect in the coupling regime relevant 
to most CDW materials. 

In this paper we extend the variational scheme proposed for the spinless Holstein model 
at half-filling by H. Zheng, D. Feinberg and M. Avignon 32 (ZFA) to finite temperature. 
We note that at zero temperature the ZFA results are in good agreement with different 
numerical works. 33-35 Actually this approach introduces lattice fluctuations on the Peierls 
dimerization since it takes into account the nonadiabatic polaron formation. Within this 
scheme all electrons in the Fermi sea are involved in the el—ph scattering leading, above the 
CO transition, to a stable phase of disordered small polarons (SP) in the intermediate to 
strong coupling regime. Thus, at half-filling the phase diagram is consistent with previous 
DMFT results pointing out a strong dependence of the transition temperature on the 
adiabaticity particularly in the weak to intermediate coupling regime. 29-31 Within the ZFA 
scheme we have derived the spectral properties useful to characterize the different phases of 
the system. Next, by using the formalism of generalized Matsubara Green's functions 36-40 , 
we have determined the scattering rate of the quasi-particles finding that, in the CO phase, 
the single-phonon emission and absorption represent the main mechanism of damping. The 
lattice fluctuation effects included in the scattering rate are beyond the ZFA scheme and 
modify the density of states that, this way, is able to capture the features of the intermediate 
coupling regime in agreement with previous DMFT studies. 29-31 

The scattering rate turns out fundamental also to derive the optical properties of the 
system. With decreasing T or increasing the el — ph coupling, our conductivity spectra are 
characterized by a transfer of spectral weight from low to high energies with the opening of 
an optical gap. The effect of the lattice fluctuations is able to smooth the inverse square-root 
singularity of the mean-field approach in the optical absorption and induces a subgap tail. 
In the weak to intermediate coupling regime, a two-peak structure is obtained within the 
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CO phase at low frequency: the remnant Drude-like term and the incipient absorption band 
centered around the gap energy. 

In section I the model and the variational approach are reviewed; in section II the spectral 
properties are discussed; in section III the damping of the particle motion is calculated; in 
section IV the optical properties are deduced. 



I. MODEL AND VARIATIONAL APPROACH 

In this paper we study the spinless Holstein model at half-filling. 23 The Hamiltonian 

reads 



H = -t g cjcj + UJ a l a i + 9^0 C \°i ( a i + a i) 
<i,j> i i 

-/iYJcfci. (1) 



Here t is the electron transfer integral between nearest neighbor (nn) sites < i,j >, cj (q) 
creates (destroys) an electron at the i-th site and \x is the chemical potential. In the second 
term of eq.(l) a] (a*) is the creation (annihilation) phonon operator at the site i, u = ^Jk/M 
denotes the frequency of the optical phonon mode, with k restoring force per length unit and 
M mass of the local oscillator. The parameter g represents the coupling constant between 
electrons and phonons 

9 = fol\/f ? ( 2 ) 

where A is the energy per displacement unit due to the coupling of the charge carriers with 
the lattice. The dimensionless parameter A 

A = ^, (3) 

with D dimensionality of the system, indicating the ratio between the SP binding energy 
and the energy gain of an itinerant electron on a rigid lattice, is useful to measure the 
strength of the el — ph interaction in the adiabatic regime. We consider spinless electrons 



since they, even if at a very rough level, mimic the action of an on-site Coulomb repulsion 
preventing the formation of local pairs. 

The hopping of electrons is supposed to take place between the equivalent nn sites of 
a simple cubic lattice separated by the distance \n — n'\ = a. The units are such that the 
Planck constant h — 1, the Boltzmann constant /c_b=1 and the lattice parameter a—1. 

Following the ZFA variational scheme, 32 we perform three successive canonical transfor- 
mations to treat the electron-phonon interaction variationally and to introduce the charge- 
ordering solution. The first is the variational Lang-Firsov unitary transformation 32 



U\ = exp 



9 (f4 c 3 + A i) ( a i _ 



(4) 



where / and Aj are variational parameters. The quantity / controls the degree of the 
polaron effect and Aj denotes a displacement field describing lattice distortions due to the 
average electron motion. The second transformation is 



U 2 = exp 



O 



(5) 



where the variational parameter a determines a phonon frequency renormalization. The 
transformed Hamiltonian H = U2 1 U{ 1 HUiU 2 is 

H = -t X i X j C Uj + ^0 E a i a i + LUJ Sillh2 ( 2C + 9 2<jJ S A< i 

<i,j> i i 

+uj sinh (2a) cosh (2a) ^ (a\a\ + aia^j — guj e 2a ^ Aj (a,i + afj 

i i 

+gu (1 - /) e 2a Y c\d (at + a}) + E c Ui (r)i ~ , (6) 



where we have the phonon operator 



Xi = exp [gfe 2a (a, - a}) 



(7) 



the renormalized phonon frequency u = u cosh(4a), the number of lattice sites L and the 
quantity rji 



Vi = <?W {f - 2) + 2^ 2 cu (/ - 1) Aj 



(8) 



At half-filling the charge-ordered solution is obtained by assuming 

A i = A + Acoe M , (9) 
where A represents the lattice distortion unaffected by the instantaneous position of electrons 

— * 

and Aco the additional local lattice distortion due to the Peierls dimerization with Q = 

(7T,7T,7r). 

In the ZFA approach the free energy is deduced employing the Bogoliubov inequality 
and introducing a test Hamiltonian characterized by non interacting electron and phonon 
degrees of freedom such that (H —H test ) t = 0, where <> t indicates a thermodynamic average 
obtained by using H test . The test Hamiltonian is given by 

Htest = -t e ff E c l c j + E ahi + Luj sinh 2 (2a) + Lg 2 uj (A 2 + A 2 co j 

<i,j> i 

-2g 2 LU A co (1 - /) E e&*4a - ^ E ^ ( 10 ) 

i i 

where the subsidiary chemical potential /j,q is 

fio = H- g 2 oj f (/ - 2) + 2g 2 co (f - 1) A. (11) 

The quantity t e ff = te~ Sr denotes the effective transfer integral, where the quantity 

St = 9 2 f 2 e~ 4a (2N + 1) (12) 

controls the band renormalization due to the antiadiabatic polaron effect, with N the aver- 
age number of phonons with frequency ujq. The test Hamiltonian is diagonalized by a third 
canonical Bogoliubov transformation 32 yealding 

Htest = E (£k +) - fM>) 4 rf k + E (f k _) - vo) pIp* + E a i a i + 

ktNZ keNZ q 

+Lu sinh 2 (2a) + Lg 2 u (A 2 + A 2 co ) , (13) 

where c4 (d k ) creates (destroys) a quasi-particle in the upper band £ k + ^ = \J^ + E 2 , p\. (p k ) 
creates (destroys) a quasi-particle in the lower band £ k ^ = —\Jt\ + E 2 and e k is the polaronic 
band. We note that NZ indicates the New Brillouin Zone defined by the condition e k < 0. 
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In the charge-ordered phase a gap opens between the upper and lower bands in the quasi- 
particle spectrum and it is twice the quantity 

E = 2g 2 u A co (!-/)■ (14) 

At half-filling (p = 0.5), the variational free energy of the system is 

j = f el + T log (l - e~^°) + co sinh 2 (2a) + g 2 io (A 2 + A 2 co ) + p(p)p (15) 

where f3 is the inverse of the temperature. The electron free energy reads 

f el = T [°. deg(e) log \n F (^)n F (^)} , (16) 
J —w L J 

where g(e) is the density of states, W = 6t e // is the renormalized band half-width, np{e) 
is the Fermi distribution, = Ve 2 + E 2 and = — ^ ( - + - ) . In order to simplify the 



calculations, in the three dimensional case we consider a semicircular density of states 



9(e) 



W 2 



0(W -\e\)\/w 2 -e 2 , 



(17) 



where 9{x) is the Heaviside function. Actually g (e) represents a simple approximate expres- 
sion for the exact density of states and it is generally used for a 3-D lattice. 42 ' 43 
Within the variational approach the kinetic energy mean value < T > 

J2 



< T >-- 



_ deg(e)- 



l-w yv ' Ve 2 + E 2 
and the electron order parameter m e 



nF (£(+))_ nj? (£i-;) 



(-)^ 



(18) 



E I" de- ^ 



(+)^ 



(19) 



D V J-w x'< 2 ■ E 2 

can be evaluated. The CO phase is characterized by the order parameter different from 



zero. 



In this paper we discuss results valid in three dimensions. In Fig.l we report the phase 
diagram obtained within our approach for W = 6ujq, with W the bare band half-width. We 
compare it with the mean-field result obtained in the adiabatic limit for / = and a = 



that overestimates the transition temperature. The A phase represents the Fermi-liquid- 
like normal state, while the B phase the disordered localized SP normal state. This latter 
phase is determined when there is absence of order (Aco = 0) and the lattice presents the 
largest distortions (/ = 1). While the transition from A to CO is continuous, the transition 
from CO to B and the crossover from A to B is rather discontinuous. We stress that a 
discontinuous character for the transition can be considered as a drawback of the ZFA 
approach. 32 However, the phase diagram bears a strong resemblance to that derived with 
more sophisticated tecniques. 29 ' 30 

In the inset of Fig.l, the ratio of T = gap A to the ordering temperature Tco is 
shown. It is interesting to note that within the CO phase it is possible to distinguish two 
regimes. The first is the mean-field regime in the sense that the ratio assumes values near 
the BCS A/ksTco = 3.52, while the second regime exhibits larger values of this quantity. 
For instance at A = 0.9 the ratio A/ksTco is about 13, the value measured in nickelates. 12 
Hence, a strong coupling CO phase develops starting from intermediate values of the el — ph 
coupling. This CO regime is characterized by a gap of the order of 2g 2 u and by distortions 
comparable to those of the disordered SP phase. 

Previous DM FT results 31 have pointed out a strong dependence of the transition temper- 
ature on the adiabaticity particularly in the weak to intermediate coupling regime. Therefore 
we have determined the transition temperature for two different values of the adiabaticity 
ratio (see Fig. 2a ). The temperature T C o strongly depends on the adiabaticity for interme- 
diate el — ph couplings. Here the lattice fluctuations are very effective in reducing Tco- The 
variation ATco = T co — Tco of the transition temperatures due to different masses of the 
oscillators is reported in Fig. 2b. The ratio between the old and the new masses simulates the 
isotopic substitution of the oxygen from 16 (M) to 18 (M*). Actually in many perovskite 
oxides the presence of the oxygen atoms induces the el — ph coupling. 44 Therefore we can 



consider a change of the values of u and g to ujq = ujq^J M/M* and g* = g (M* /M) 1 ^ 4 
respectively. 45 The softening of the phonon frequency induces an increase of the transition 
temperature. Furthemore, as a function of A, the difference between the two ordering tern- 



peratures is characterized by a maximum in the region of weak to intermediate el — ph 
couplings. 



II. SPECTRAL PROPERTIES 



In this section we calculate the spectral properties within the ZFA approach. They are 
discussed in order to characterize the different phases of the diagram reported in Fig.l. 

The electron Matsubara Green's function can be disentangled into electronic and 
phononic terms 46 by using the test Hamiltonian (10), hence 



G (k, r) = - (i) £ e ik -( Ri - R ^) (T t q (r) c)) t (T T N (r) X}) u 



(20) 



where we have Q (r) = e rHt Qe rHt and Xi(r) = e THt Xie rHt . The Green's function in 
Matsubara frequencies u n becomes 

G(k,iu n ) = e- ST G [c) (k,ioj n ) + 

e ~ ST (^) E (^) E^ W (k u tuj m ) m )-{ e *cosh[a (r-f)] _ 1} (21) 

where s is 



s = 2/^ 2 e- 4a [N (N + 1)] . 



(22) 



In eq.(21) (k, iu n ) is the Green's function 



S (c) (k, i Wn ) = (k, iu n ) + v j£H (k, mo n ) , 



(23) 



where (k, iuj n ) and C/( ) (k, iu n ) are the Green's functions associated to the quasi- 
particles of the upper and lower bands, with -u 2 . given by 



1 + 



ek 



and by 
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(24) 



(25) 



Two physically distinct terms 46 appear in the Eq.(21). The first derives from the coherent 
motion of quasi-particles and their surrounding phonon cloud, while the incoherent term is 
due to processes changing the number of phonons during the hopping of the charges. 

Making the analytic continuation iu n — > u+iS, one obtains the retarted Green's function 
G ret (k, uj) and the spectral function 



A(k, u) = -2ZG ret (k, oo) = 2ixe~ s ? [u 2 k 5 (u - + v 2 5 (u - & +) ) 



+ 



+°° m- 



2ne- ST h(s)e—[l - n F (u - lu) )}K(uo - lQ Q ) + 

Z(^0)-oo 

2ixe~ ST Ii(s)e~^n F (uJ + lO )K(uj + lQ ) + 

Z(^0)-oo 

2ne~ ST [I (s)-l]K(Lu), (26) 
where the function K{yj) is 

K(u) = H(u) + H(-cu), (27) 



with 



H{uo) = V e2 1 Q (uj-E)e (VE 2 + W 2 



E , v -•)• (28) 
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In the normal phase the first term of eq.(26) represents the purely polaronic band contribu- 
tion and shows a delta behavior. In the CO phase, this term is equal to the result of the 
mean-field approach when one neglects the renormalization of the upper and lower band due 
to the polaron effect. The remaining terms of eq.(26) provide the incoherent contribution 
that spreads over a wide energy range. 

Thus we get the renormalized density of states N{uj) that is normalized to unity. In 
Fig. 3a we report the density of states for three different values of the el — ph coupling at 
T = 0.023VT. It is apparent that in the normal phase (A = 0.15) the weight of the coherent 
term is prevalent and the density of states is practically unchanged with respect to the non- 
interacting case. Entering the CO phase, a gap opens in the quasi-particle spectrum and it 
broadens with increasing the parameter A. We note that, at the energies corresponding to 
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the gap, the inverse square-root singularity occurs. The other sharp features of the density 
of states are due to one-phonon and two-phonon processes in the upper and lower bands. 

At higher temperatures (Fig. 3b) similar features are found. However, in this case, for 
large el — ph couplings, the density of states is determined by the incoherent dynamics of 
the small polaron excitations. 46 It is made up of two bands peaked approximatively around 
—g 2 iVo and g 2 u>o, whose heights are equal at half-filling. We stress the evolution of the gap 
energy of the CO phase toward the characteristic energy 2g 2 u that represents the energy 
range between the two peaks of the SP density of states. 



III. DAMPING 

In this section we deal with the el —ph self-energy that includes lattice fluctuation effects 
beyond the ZFA approach. This quantity allows to determine the scattering rate of the 
quasi-particles of the upper and lower bands and a better density of states. This analysis 
will play an essential role in the infrared absorption calculations that is the main aim of this 
paper. For sake of clarity we will put forward this discussion. 

Retaining only the dominant autocorrelation terms at the second step of iteration 36-40 , 
we derive the local self-energy 

S (2) (iu; B ) = Q E (^J E£ (C) (ki,^n) [ ' dre^-^iMr) + / 2 (r)], (29) 



where /i(r) is 



A(r) = Zt 2 e 



2 -2S 7 



3 2cosh[o;o(T"-|)] 



(30) 



z = 2s, with s given by eq.(22), Z the number of nn sites, and /2(t) is 



/ 2 (r) = 2^ V a (l - f) 2 [N (N + lp cosh 



UJ T 



S 



(31) 



Other terms beyond the autocorrelation ones induce a dependence of the self-energy on the 
momentum and could be included. 40 We note that the fundamental autocorrelation contribu- 
tion presents some analogies with the corresponding self-energy within DM FT approaches 
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where the independence on the momentum is the main assumption when the interaction 
term is local. 

Making the analytic continuation iuj n — > u + i5, we can calculate the scattering rate of 
the quasi-particles of the upper band 



r + (k) = f( w = t y ) = -29e£> ( w = • 



We can express the scattering rate in the following way 



T+(k) — r + (£fc + ^) — Ti^ phon (^') + T multi-phonic ') + FresiC), 

where Ti_ phon is the contribution due to single phonon processes only 

ri- pfto „(d +) ) = 2Zt'e-^h(z) sinh 9l , l=1 (& ] ) + 

Fmuiti-phon represents the scattering rate by multiphonon processes 



*(+) 



(+)^ 



1 multi—phon 



/to f 



= 2Z*e-** £ (*) -h ( ^ ) ^(£ k +) ) 



(32) 



(33) 



(34) 



(35) 



and T res denotes a residue term due to the difference between the coherent and the incoherent 
contribution 47 



r P «(d +) ) = Zt 2 e~ 2S - [Jo (*) - 1] B(& } ). 



(+)^ 



In the previous equations the function 5 f i,z(Ck + ' > ) rea ds 



<7i,/(£k +) ) = ^o(to ) + n^' + too) 5(4 +J + too) 



(+) 



+ 



N (lu ) + 1 - n F (£ k +) - too) B(&> ~ to ) 



and # 2 (&" 



iV + n F (£ k +) +^ ) 5(^+^o) + iV + l-n F (^-^o) 5(^-^o) 



:(+) 



(+) 



(36) 



(37) 



(38) 
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Moreover the function B(x) is equal to 2irK(x), with K(x) given in eq.(27). The decompo- 
sition of the scattering rate in three distinct terms has been introduced in order to simplify 
the analysis of the results. 

We define the scattering rate for the quasi-particles of the lower band 

r_(k) = f ( W = = -29E£> ( w = , (39) 

that turns out to be equal to T + (k). Thus we can take into account only one scattering rate 

r(k) = r + (k) = r_(k). (40) 

In Fig. 4a we report the scattering rate T for different el—ph couplings at low temperature. 
In the normal phase (A = 0.15) the scattering rate is zero within u of the chemical potential 
\i. Therefore at this temperature the main mechanism of energy loss is the spontaneous 
emission of phonons with frequency u by the quasi-particles (one phonon processes are 
prevalent). Thus the behavior of the scattering rate is determined by of the Fermi statistics 
since the quasi-particle excitations within u of the Fermi energy cannot lose energy 46,47 In 
the CO phase (A = 0.50) the gap in the scattering rate opens at energies given by the sum 
of the gap energy (indicated by the arrows in the figure) and the phonon frequency. Since 
there are not available states within the gap, only the quasi-particles with such energies can 
be scattered by the residue el — ph interaction. 

In Fig. 4b we concentrate on the scattering rate at a fixed el — ph coupling (A = 0.4) 
for different temperatures. At low T (T = 0.01VF) there is the gap due to the Peierls 
dimerization and the processes of phonon spontaneous emission by the quasi-particles. With 
rising temperature, we can follow the increase of the quantity T in the normal phase because 
of the absorption and emission of phonons. In the normal phase at finite temperatures, these 
scattering processes are effective also in the energy range around the chemical potential where 
there are able to enhance the rate. Actually, at higher temperatures (T = 0.14VF), the rate 
is weakly dependent on the energies involved in the phonon scattering. It is confirmed that 
in this regime of temperatures the single-phonon emission and absorption give an important 
contribution to the quantity T. 
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For large el — ph couplings and for small polaron excitations, the quantity T decreases 
when T increases but it is always larger than t e ff, therefore the electronic states lose their 
individual characteristics and the electron motion is predominantly a diffusive process. A 
high-temperature expansion 37 ' 38 ' 46 provides the scattering rate 



(41) 



with z = 2g 2 uj and E a = z/A. This value coincides with the well-known rate of the Holstein 
polaron in the "classical" limit, with E a typical activation energy for hopping. 23 

The introduction of the damping allows to improve the approximations of calculation 
for the spectral properties. This can be carried out substituting in eq.(21) the new Green's 
function 

(k, iw n ) = ulQ^ (k, iu n ) + vlG^ (k, iu n ) . (42) 
for Q( c \ The Green's function (k, iu n ) is expressed as 

s^ n)= r^*^}, ,43) 

where the spectral function is assumed to be 



AM(k, w ) = ^ n-, 

[r(k)]V4 + (c-^ ) ) 2 



(44) 



with v standing for + or — . We can determine a new electron spectral function and density 
of states Npi uct that includes fluctuation effects beyond the ZFA approach. These effects 
are able to change the density of states in the Fermi-liquid-like phase. Indeed in this phase 
it is interesting to study the behavior of the density as function of the temperature at a fixed 
value of the coupling constant A (see Fig. 5a). While the density obtained within the ZFA 
approach is similar to the bare one (solid line in figure), now N F i uct shows an enhancement 
at the Fermi energy in the low temperature regime (T = 0.07VT) and a depression at 
the same energy in the high T case (T = 0.20VT). Actually at low T the processes of 
phonon spontaneous emission subtract spectral weight to states of energy larger than c^o 
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compared with the chemical potential inducing an increase at the Fermi energies where the 
states are not damped. At intermediate temperatures (T = 0.12W), where real phonons 
are present, the absorption and emission processes are effective also around the chemical 
potential and oppose this enhancement, so that the density of states at the Fermi energy is 
nearly unchanged. At higher temperatures these processes dominate causing a reduction of 
the spectral weight at the chemical potential. 

In Fig. 5b the density of states N F i uct is plotted as a function of the coupling constant A 
showing that it is able to capture the features of the intermediate coupling regime. For these 
values of parameters, the density within the ZFA approach shows a discontinuous transition 
from a lightly renormalized function to a strongly modified density with a marked reduction 
at the chemical potential. Due to the included lattice fluctuations, N F i uct evolves gradually 
from the bare density to the characteristic small polaron function showing pseudogap features 
in agreement with previous studies. 29-31 Not only in the normal phase but also in the CO 
state, the quantity Y changes the spectral properties. Actually, in the boundary phase region 
below the transition temperature of the ZFA approach, a pseudogap opens in the density 
of states as a precursor effect of a gap at lower temperatures and stronger el — ph couplings. 

Thus the lattice fluctuation effects introduced by means a self-energy insertion turn out 
to be able to correct the drawbacks of the ZFA approach. Next we will focus on the optical 
properties analyzing the prominent role of the lattice fluctuations in determining the infrared 
absorption of the system. 



In this section we deal with the optical properties near and within the CO phase. Cal- 
culating the real part of the conductivity in a regime of linear response reduces to evaluate 
the retarded current-current correlation function 



IV. OPTICAL PROPERTIES 




(45) 
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The electron motion will be supposed to take place between the equivalent nn sites of the 
cubic lattice, hence the tensor a a ^ is assumed to be diagonal with mutually equal elements 



' a,a- 



Performing the two canonical transformations (4,5) and making the decoupling 46 of the 
correlation function in the electron and phonon terms through the introduction of H test (10), 
we get in Matsubara frequencies 

n Q , Q (^ n ) = e 2 t 2 (-1) £ £(* • 5') f dre**> T $(i, i f , 5, 5', r)A(i, i', 5, 5', r), (46) 

where the function A denotes the electron correlation function 



A (i, i', 5, 5', r) = (T T c\(T)c i+ s & (T)4' + s> & c i >) t 



and the function $ the phonon correlation function 



(47) 



$ (i, i', 5, 5>, r) = (T T A > /(r)A > t+ ^(r)4 + ,, (i X i ,)i. 



(48) 



To derive the optical properties, the role of the damping T of the particle motion is 
fundamental. Since the electron correlation function can be expressed as a function of the 
Green's functions the effect of the damping 46 can enter our calculation substituting 
for given by eq.(43). Furthermore, in order to simplify the analysis of our results, 
we separate $ into two contributions 



W) = [(X}x i+S& ) t ] 2 {<f>(i,i',6,6') - [(XjX l+s& ) t ] 2 } 
^(t,t',5,5')-e- 2ST . 



e -2S T + 



(49) 



Considering the two terms of eq.(49), the current- current correlation function can be written 

as 



n a ^(<a;0 = nW(ia; B ) + n«(<a; B ). 



(50) 



The first term reads 



n«(^ n ) = 4e 2 tV 2 ^ (j) £ sin 2 (A; Q ) £ (k) S<^>(k,ia; B ), (51) 
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where we have 



(k) =/>(->-) (k) = (*4-^) 2 , (52) 

h i+ ^ (k) = (k) = 4u£u£ (53) 

with S'(^>^)(k,iw n ) given by 

S {ui ' U2 \k,tuj n ) = [ P dre tuJnT Q (ui \k,-T)g^\k,T). (54) 
Jo 

The second term of the current-current correlation function is obtained retaining only the 
main autocorrelation term % = %' and 5 = 5' 

ngi(^n) = (^f) (jj) E E EP M (Mi)^ >(Mi,™»), (55) 
where we have 

(k, kx) = (k, kO = ( Mk « kl - t; k t; kl ) 2 , (56) 

p< + >") (k, kO = (k, kx) = (« kVkl + Mkl t; k ) 2 , (57) 

the function T^'^k, k 1; iw n ) 

T^^)(k,k 1 , ? ^ n ) = / /3 dre iw » T ^)(k,-r)^)(k 1 ,r)/ 1 (r), (58) 

with /i(r) given by eq.(30). 

We perform the analytic continuation iu n ^ u + i5, and, clearly, the conductivity can 
be expressed as a sum of two terms 46 



9 



$la a ,a M = L — J = 9fo7& h ) (w) + $la^ (u). (59) 



a; 



As in the spectral properties, the appearance of two physically distinct contributions, the 
coherent and incoherent one, occurs. Actually the first term tfta^J^ is due to the charge 
transfer affected by the interactions with the lattice but not accompanied by processes 

17 



changing the number of phonons. On the other hand, the incoherent term 9?<7^™ coh ) in eq. 
(59) derives from inelastic scattering processes of emission and absorption of phonons. The 
coherent conductivity is derived as 

'4e 2 t 2 \ ro 



e" 25T E / rfe[n F (e (iyi) -cu)-n F (e (iyi) )]C ,(l/1 ' iy2) (e,^)/i(e)^ 1 '" 2) (e), 

"1,1*2 J ~ W 



where C^ U2 \e,u) is 



Tie) 



T 2 {e) + {^ +uj f 



and h(e) is defined as 



In eq.(60) the function A^ 1,U2 \e) is expressed by 



(60) 



(61) 



(62) 



A (+ ' +) (e) = A^-\e) 



e 2 + E 2 



(63) 



and 



A<+->(e) = A(--+)(e) = 



£ 2 



e 2 + E 2 



(64) 



The latter term of the conductivity becomes 



2+2" 



2eH 



00 



-2S 7 



E [°- rfe /°- dei^(e)^(ei)i2 ( ^^ ) (e>ei^) + 

A, „n J-W J-W 



(—) e- 2ST Uz) - 1] E >< 



x /°_ rfe /°_ de^e)^) [n F (^ 2) -c) -n F (^ 2) )l C^'^(e, Cl , «), (65) 



where g(e) is the density of states (17), the function R^ ul ' U2 \e, ei,u) is given by 



'Pu l 



R^(e,e 1 ,u) = 2£/ t (*)sinh ( ^ ) J^»( £|£l)(1 ,) + flj^( £)£l)W ) 



(66) 



C^'^e^x) is 
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c*«W) = ; TO + rfa)] 

2 [r ( £ ) + r( £l )]V4 + (£<-> + 

and 



^ = v jy /e? + E*. (68) 
We notice that the functions j\ Ul,U2 \e,ei,uj) 

J^ 2) (e, e u u) = C^\e, e u oo + lcu ) [n F {^ ] - loo - u) - n F {^ - x 



x 



'N (lu )+n F (& 2) )} (69) 
and H\ v ^ 2) '(e,ei,w) 

H^\e,e 1 ,cu) = C ( ^\e, e u u - l^)[n F {^ + to - w ) - n F (^ 2) + ia> )] x 



x 



7V (to ) + l-n F (^ 2) ) (70) 



describe phonon absorption and emission processes, respectively. 

In the limit of high temperatures and for small polaron excitations, the phononic inco- 
herent absorption is prevalent and in this limit an analytic expansion can be performed. 46 

The internal consistency of our approach is verified by using the sum rule 

jf du&a ata (u) = --e 2 < T a>a > (71) 

where (T Q Q ,) is the mean value of component of the kinetic energy equal to one third of 
eq.(18). Employing the calculated spectra, we have checked tha the two sides of eq.(71) 
differ in a few per cent. 

In Fig. 6 we report the conductivity spectra for T = 0.023M 7 at el — ph couplings. With 
rising the el — ph coupling, a tranfer of spectral weight from low to high energies takes 
place. The Drude term makes smaller and, in the weak to intermediate CO regime, the 
optical response is characterized by two components: the remnant Drude-like term and the 
incipient absorption band centered around the gap energy. For stronger el — ph couplings, 
the band is peaked around the gap but it is accompained by a subgap tail due to the lattice 
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fluctuations. 19,20,48 Thus, by increasing the parameter A and by crossing the phase transition, 
we do not have any sharp changes in the optical spectra but a rather continuous evolution. 
We point out that in the mean-field approach the infrared absorption occurs only at energies 
above the gap and it is characterized by the inverse square-root singularity. 49 

In the inset of Fig. 6 the gradual disappearance of the Drude term and the rise of the 
interband absorption can be followed in a narrow range of energies. At higher temperatures 
a similar evolution of the bands of the infrared absorption is found by increasing the el — ph 
strength. Actually, in both the ranges of temperature, the response in the limit u — > is 
completely suppressed only for strong el — ph couplings where a well-defined optical gap 
opens. Furthermore in the CO phase there is a lowering of the height of the interband terms 
that spread over a wide range of energies. 

It is interesting to analyze the conductivity spectra at a fixed el—ph coupling for different 
temperatures. We take into account two different regimes: A = 0.4 (A — > CO phase) and 
A = 0.8 (B — > CO phase). In the first case (Fig. 7), by starting from the normal phase and 
by lowering the temperature, there is a redistribution of the spectral weight among the two 
components accompanied by a progressive narrowing of the Drude term. The maximum of 
the second band slightly shifts towards higher frequencies and the optical gap fully develops 
at low temperatures well within the CO phase. We stress that these results of the optical 
absorption are consistent with experimental measurements in dichalcogenides. 18 

In the other regime (A = 0.8) very different features of the optical response are obtained 
(Fig. 8). Here the evolution from the SP to CO absorption band is characterized by near 
peak energy. In the CO phase, below Tco, there is a residue of optical response in the 
limit uj — > and the optical gap starts to appear for lower temperatures in corrispondence 
of energies that can decrease with lowering T. These results can be closely connected with 
the conductivity spectra measured in nickelates 12 ' 15 that reveal the relevance of the lattice 
degrees in stabilizing the CO phase. 
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V. SUMMARY AND CONCLUSIONS 



We have discussed the phase diagram, the spectral and optical properties of the half- 
filled spinless Holstein model in three dimensions as a function of the temperature and the 
el — ph coupling. The lattice fluctuation effects play a crucial role in determining the density 
of states and the optical response of the system. 

The phase diagram is consistent with previous DMFT results pointing out a strong 
dependence of the transition temperature on the adiabaticity particularly in the weak to 
intermediate coupling regime. The density of states captures the features of the intermediate 
coupling regime in agreement with DMFT studies. 

Concerning the optical absorption, we have observed that, with decreasing T or increasing 
the el—ph coupling, the ordered phase affects the conductivity spectra inducing a transfer of 
spectral weight from low to high energies with the opening of an optical gap. The inclusion of 
effects of quantum lattice fluctuations is able to smooth the inverse square-root singularity of 
the mean-field approach and induces a subgap tail absorption. In the weak to intermediate 
coupling regime, a two-peak structure is obtained within the CO phase at low frequency: the 
remnant Drude-like term and the incipient absorption band centered around the gap energy. 
The results obtained in the intermediate and in the strong coupling regime are consistent 
with experimental conductivity spectra. 12,15,18,20 

In this paper lattice fluctuation effects beyond the variational ZFA approach are included 
by means of self-energy insertions. At the lowest order the main contribution to the self- 
energy is independent on the momentum, however at the same order other terms could be 
easily included in order to give a non local quantity. Furthermore, the scattering rate can be 
calculated not only perturbatively but also employing a self-consistent procedure. 36-38,46 If 
also the real part of the self-energy is introduced by using the Kramers-Kronig relation, the 
procedure of calculation could become more complete. The agreement of the results obtained 
in this paper with other approaches seems to suggest that a second-order perturbation theory 
on the good ZFA solution can capture the physics of the model in the different regimes of 
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the normal and ordered state. 

Finally we note that our approach is valid in the infrared range of frequencies where the 
interband absorption takes place. Thus it is not able to reveal the structures attribuited to 
collective excitation modes arising from the CDW condensate. 1,16 ' 17 ' 5 ° 

FIGURE CAPTIONS 

Fl The phase diagram (solid line) corresponding to W = 6u . CO characterizes the charge 
ordering phase, the A phase represents the Fermi-liquid-like normal state, the B phase 
the disordered small polaron normal state (the crossover between the A and B phase 
is marked by the dashed line). In the inset, the ratio of T = gap A to the ordering 
temperature Tqo as a function of A. 

F2 (a) The transition temperature as a function of A for two different values of the adia- 
baticity ratio: W = Qu (solid line) and W = 3u (dashed line), (b) The variation 
ATco = Too ~ ^co of the transition temperatures due to different masses of the 
oscillator as a function of A. 

F3 (a) The renormalized density of states at W = 6uo and T = 0.023VF as a function of 
the energy (in units of uq) for different el — ph couplings: A = 0.5 (solid line), A = 0.35 
(dashed line) and A = 0.15 (dotted line). 

(b)The renormalized density of states at W = 6u and T = 0.14VF as a function of 
the energy (in units of u ) for different el — ph couplings: A = 0.65 (solid line), A = 0.6 
(dashed line), A = 0.56 (dotted line) and A = 0.4 (dash-dotted line). 

F4 (a) The scattering rate at W = Quj and T = 0.0231F as a function of the energy (in 
units of uo) for different el — ph couplings: A = 0.5 (solid line) and A = 0.15 (dashed 
line). 

(b) The scattering rate at W = 6ujq and A = 0.4VF as a function of the energy (in 
units of uj ) for different temperatures : T = 0.01VF (solid line), T = 0.05PF (dashed 
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line), and T = O.UW (dotted line). 
The arrows indicate the gap energy. 

F5 (a) The new renormalized density of states at W = 6u and A = 0.2 as a function of 
the energy (in units of u>o) for different temperatures. 

(b)The new renormalized density of states at W = Glu and T = 0.2W as a function 
of the energy (in units of c^o) for different el — ph couplings. 

F6 The conductivity up to 6 u at different el — ph couplings. In the inset the conductivity 
up to 2 ujq at different el — ph couplings. The conductivities are expressed in units of 
, with m = 757. 

F7 The conductivity (in units of with m = ^) up to 2 lo at different el—ph couplings. 
F8 The conductivity (in units of with m = ^) up to 16 uj at different el—ph couplings. 
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